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ABSTRACT 

We report on X-ray spectroscopic observations with XMM/Newton of the ul- 
tracompact, double white dwarf binary, GP Com. With the Reflection Grating 
Spectrometers (RGS) we detect the La and L/3 lines of hydrogen-like nitrogen 
(N VII) and neon (Ne X), as well as the helium-like triplets (N VI and Ne IX) 
of these same elements. All the emission lines are unresolved. These are the first 
detections of X-ray emission lines from a double-degenerate, AM CVn system. 
We detect the resonance (r) and intercombination (i) lines of the N VI triplet, but 
not the forbidden (f) line. The implied line ratios for N VI, R = f ji < 0.3, and 
£r=(/ + z)/r«l, combined with the strong resonance line are consistent with 
formation in a dense, collision-dominated plasma. Both the RGS and EPIC/MOS 
spectra are well fit by emission from an optically thin thermal plasma with an 
emission measure (EM) oc (kT/6.5 keV) as (model cevmkl in XSPEC). Helium, 
nitrogen, oxygen and neon are required to adequately model the spectrum, how- 
ever, the inclusion of sulphur and iron further improves the fit, suggesting these 
elements may also be present at low abundance. We confirm in the X-rays the un- 
derabundance of both carbon and oxygen relative to nitrogen, first deduced from 
optical spectroscopy by Marsh et al. The average X-ray luminosity of ~ 3 x 10 30 
ergs s -1 implies a mass accretion rate m ~ 9 x 10~ 13 M Q yr _1 . The implied 
temperature and density of the emitting plasma, combined with the presence of 
narrow emission lines and the low m value, are consistent with production of the 
X-ray emission in an optically thin boundary layer just above the surface of the 
white dwarf. 



Subject headings: Binaries: general - Stars: individual (GP Com) - Stars: white 
dwarfs - cataclysmic variables - X-rays: stars - X-rays: binaries 
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1. Introduction 

The AM CVn stars are among the most compact binary systems known. Their optical 
spectra are dominated by broad helium emission lines originating in an accretion disk formed 
by mass transfer from a degenerate helium dwarf onto the primary white dwarf. Their 
orbital periods range from about 10 minutes to 1 hour. They are natural laboratories for the 
study of poorly understood binary evolution processes, such as common envelope evolution. 
Moreover, the absence of hydrogen allows for the study of accretion disks dominated by 
helium. They may also be a significant channel for the production of Type la supernovae 
and neutron stars via accretion induced collapse (for a review see Warner 1995). 

GP Com is one of the better studied, nearby AM CVn systems. It has a 46.5 minute or- 
bital period (Nather Robinson & Stover 1981), resides at a distance of i=» 70 pc (Thorstensen 
2003), and has been studied extensively in the optical and UV. Recent studies by Marsh 
(1999) and Morales- Rueda et al.(2003) have used optical and UV+optical spectroscopy to 
probe the dynamics of the system in great detail. In addition to helium, nitrogen emission 
lines have been seen in the UV (Lambert & Slovak 1981; Marsh et al. 1995), and optical 
(Marsh, Home & Rosen 1991). These observations identified a nitrogen overabundance, rela- 
tive to carbon and oxygen, as well as an underabundance of heavy elements (for example, Ca, 
Si and Fe). The low metallicity may be consistent with its suggested halo origin (see Giclas, 
Burnham & Thomas 1961). The evolution scenario favored by Marsh et al. (1991) posits an 
initially low metallicity for the progenitor stars. The CNO elements are subsequently pro- 
duced within the primary star and then mixed throughout the secondary during a common 
envelope phase. The nitrogen abundance is then increased by CNO-cycle hydrogen burning 
within the secondary (ie. the mass donor). Finally, a number of ultracompact systems with 
neutron star primaries also show apparent neon enrichment (see Juett & Chakrabarty 2003; 
Schulz et al. 2001). 

Accretion should make such systems soft X-ray sources, and indeed, several have been 
detected in the X-ray band, including the prototype AM CVn, CR Boo, and GP Com (van 
Teeseling & Verbunt 1994; Ulla 1995; and Eracleous, Halpern & Patterson 1991). The X-ray 
flux from accreting, non-magnetic white dwarfs is thought to originate in a boundary layer 
that shocks and decelerates the Keplerian flow, allowing it to eventually settle on the white 
dwarf. Early work suggested that the boundary layer should be optically thick-radiating in 
the soft X-ray band and extreme UV-at mass accretion rates above mwl.6x 1O~ 1O M yr _1 
(Pringle & Savonije 1979). At lower rates the boundary layer should become more radially 
extended and optically thin, producing a harder X-ray spectrum. Aspects of this basic 
scenario have been confirmed by more recent calculations (see, for example, Patterson & 
Raymond 1985; Narayan & Popham 1993; Popham & Narayan 1995). 
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To date, X-ray observations of most AM CVn systems have been made at relatively low 
signal to noise levels and with modest spectral resolution. For example, ROSAT observations 
of GP Com found the source at a flux level ~ 1.2 x 10 -11 ergs cm -2 s _1 in the 0.2 - 2 keV 
band, and confirmed that it is variable. These observations also found that the 0.2 - 2 
keV spectrum cannot be described by a single temperature component (van Teeseling & 
Verbunt 1994), but detailed modelling was problematic due to the relatively poor statistics 
and restricted bandpass. 

In this Letter we summarize results of recent, high signal to noise X-ray spectral mea- 
surements of GP Com with XMM/Newton. These data represent perhaps the best X-ray 
spectral measurements of an AM CVn system to date, and also show the first X-ray emission 
lines from such a system. 

2. XMM Observations 

XMM/Newton observed GP Com for 56.5 ksec beginning on January 3, 2001. The 
EPIC observations were done in imaging mode. We used version 5.4.1 of the XMM software 
(SAS) to analyze the data. GP Com was easily detected at the few counts s -1 level in the 
EPIC instruments. A lightcurve combining data from the PN and MOS detectors is shown 
in Figure 1, and confirms substantial X-ray variability, including several 2 to 4 minute flares 
during which the X-ray flux increased by factors of 2 - 3. Flaring has often been seen in the 
UV and optical, particularly in the N V (1240A) resonance line (Marsh et al. 1995). We will 
present a detailed timing study of GP Com in a subsequent paper. We extracted spectra 
and produced response files for each of the instruments using the relevant SAS tools. 

2.1. Spectral Analysis: RGS and EPIC/MOS 

The spectrum extracted from both RGS instruments is shown in Figure 2, along with 
identifications for all the strongest line features. The strongest detected feature is the Lyo; 
line of hydrogenic nitrogen (N VII) at 24.78 A, but other lines of nitrogen are also present, 
including the Ly/3 line as well as the helium-like triplet of N VI at ~ 29A The only other 
signifcant line features are due to neon. The hydrogenic Lya and Ly/3 lines as well as the 
helium-like triplet are detected. Particularly striking is the absence of any strong carbon or 
oxygen lines in the spectrum. 

We fit the observed line features using Lorentzian profiles with the widths fixed at that of 
the RGS line spread function (LSF) for that wavelength (see den Herder et al. 2001; Pollock 
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et al. 2003). All the hydrogenic lines are well fit with such a profile and are unresolved. 
Table 1 summarizes the salient properties of all the detected lines. 

As is well known, the helium-like triplets can be used as probes of the plasma where the 
lines are produced (see Gabriel & Jordan 1969; Porquet & Dubau 2000). The line flux ratios, 
R = f/i, and G — (i + /)/r, of the resonance (r), intercombination (i) and forbidden (f) 
transitions are sensitive to the electron density and temperature of the plasma, respectively. 
For each triplet we fit a Lorentzian profile at the rest wavelength of each transition. To 
determine the line ratios, we fixed the wavelength and width of each line, but allowed the 
line fluxes to vary. The signal to noise ratio in both triplets is not particularly high, but 
the line ratios can at least be roughly estimated. The nitrogen triplet is better resolved 
than that of neon. The results of our fits are shown in Figure 3 (see also Table 1). For the 
N VI triplet we detect the resonance and intercombination lines, but interestingly, not the 
forbidden line. The line ratios for N VI, R = f/i < 0.3, and G = (/ + i)/r«l, and the 
strong resonance line are consistent with formation in a dense, collision-dominated plasma 
(Porquet & Dubau 2000). For nitrogen, a strong UV radiation field can also decrease the 
strength of the forbidden line, mimicking a density effect, however, the line ratios from the 
neon triplet-which are much less sensitive to radiation field effects-are also consistent with 
a plasma at or near the critical density n c = 5.3 x 10 9 cm -3 for nitrogen. 

We fit the EPIC/MOS and RGS spectra with several collisional ionization equilibrium 
(CIE) plasma emission models. We found that multi-temperature models fit the data sig- 
nificantly better than single temperature models. For example, the multi-temperature and 
variable abundance model cevmkl in XSPEC (version 11.3) provides an excellent fit to the 
data from both the RGS and EPIC instruments (see Mewe, Gronenschild & van den Oord 
1985; Mewe, Lemen & van den Oord 1986; Liedahl, Osterheld & Goldstein 1995). This 
model uses a power law in temperature for the emission measure distribution (EM), vis., 
EM oc (kT/kTmax) 1 , where kT max and 7 are the maximum temperature and power law 
index, respectively. 

There is some evidence for modest pile-up in the PN spectra, so for the broad-band 
modeling we focused on the MOS data only. Figure 4 shows our best fit to the EPIC/MOS 
data using the cevmkl model. The model fits the data well, with a minimum x 2 of 1,005 
for 1,045 degrees of freedom. Using the MOS data only, the best-fit emission measure 
distribution is characterized by kT max = 6.3 ± 0.3 keV, and 7 = 0.85 ± 0.05. The 0.3 - 
10 keV fluxes are 4.7 and 5.3 x 10 -12 ergs cm~ 2 s _1 in the MOS1 and MOS2, respectively. 
Taking the mean, we arrive at an average luminosity (at 70 pc) of 2.94 x 10 30 ergs s _1 in the 
same band. 

The elements strongly required by the fit are helium, nitrogen, oxygen and neon, how- 
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ever, x 2 additionally improves by 13.35 and 23.9 with the inclusion of iron and sulphur, 
respectively, suggesting these elements may also be present at low abundance. No other 
elements are required to fit the spectrum. Using this model we determined the abundances 
relative to the solar values tabulated by Anders & Grevesse (1989). Expressed as mass 
fractions we find; X He = 0.977 ± 0.002, X N = 1.7 ± 0.1 x 10" 2 , X = 2.2 ± 0.3 x 10~ 3 , 
X Ne = 3.7 ± 0.2 x 10~ 3 , X s = 2.3 ± 0.6 x 10~ 4 , and I Fe = 8±2x 10~ 5 . In addition, we 
placed limits (3<x) on several prominent elements which were not detected; X c < 2 x 10~ 3 , 
X Mg < 1.6 x 10~ 4 and X S i < 8 x 10~ 5 . Our limits on some of the heavier metals (such as 
Si, Mg, and Fe) suggest underabundances compared with solar values. This is qualitatively 
similar to the abundances derived from optical spectroscopy (Marsh et al. 1991). 



3. Discussion 

The nitrogen lines in the X-ray band provide additional confirmation of the overabun- 
dance of nitrogen relative to carbon and oxygen deduced from the optical (Marsh et al. 
1991). The high nitrogen abundance likely reflects the operation of CNO-cycle thermonu- 
clear processing in the secondary (Marsh et al. 1991). We find a nitrogen to helium ratio, 
by number, of ~ 4.8 x 10~ 3 , which is moderately higher than the value derived by Marsh 
et al. (1991) of 2.7 x 10~ 3 (see their Table 4). Our derived oxygen to helium number ra- 
tio is, however, about a factor of 10 higher than that deduced from the optical. Although 
systematic modelling effects could be present, these results suggest that the X-ray emitting 
gas may have higher oxygen abundance than the accretion disk (which produces the optical 
emission lines). If the accretor were an oxygen-rich white dwarf, and some of its matter is 
mixed with the accreted material in the turbulent boundary layer, then this could perhaps 
increase the oxygen abundance relative to the donor material. 

The neon lines provide the first definitive detection of an element heavier than oxygen 
in GP Com. Our derived neon mass fraction (fa 0.0037) is a factor of 2 larger than the 
solar value. Moreover, modelling of the optical spectra also suggests a neon enrichment 
of about 1.3 compared to solar abundances (Marsh et al. 1991). Neon is a by-product of 
helium burning under a wide range of conditions (see, for example, Clayton 1983; Iben & 
Truran 1983), however, the low carbon abundance is likely incompatible with helium burning 
having ocurred in the secondary. Therefore, this avenue for the neon production appears to 
be a dead-end. Again, mixing of neon from the accretor (which likely has undergone helium 
burning) into the donor material in the boundary layer could perhaps increase the neon 
abundance. It appears that gravitational settling and fractionation of neon (in the form of 
22 Ne) are also not relevant, since nitrogen is depleted in the production of this isotope, and 
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we clearly measure a large nitrogen abundance (see, for example, Bildsten & Hall 2001). 

Neon enhancements have been suggested for several ultracompact neutron star binaries; 
4U 1543-624 and 2S 0918-549 (Juett & Chakrabarty 2003), and 4U 1627-67 (Schulz et al. 
2001). For 4U 1543-624, Juett & Chakrabarty (2003) deduced Ne/O « 1.5, by number, from 
studies of absorption edges in Chandra high resolution spectra, however, we note that these 
measurements can be influenced by the unknown ionization structure in the circumbinary 
material. We find a similar value of Ne/O = 1.4 (by number) for GP Com. Yungelson, 
Nelemans & van den Heuvel (2002) suggest fractionation of neon in the white dwarf donor 
as a possible explanation for the neon enhancements in the neutron star systems. This 
scenario again requires a helium burning donor and thus seems not to work in the case of 



The X-ray measurements provide a means to estimate the mass accretion rate, m, 
onto the white dwarf primary in GP Com. Assuming matter falls from the inner Lagrange 
point onto the primary, and that half the gravitational potential energy is dissipated in the 
accretion disk, the X-ray luminosity, L x can be approximated as (see, for example, Nelemans, 
Yungelson & Portegies Zwart 2004), 



where Mi, R, R^, and m are the primary mass, primary radius, the distance from the 
inner Lagrange point to the center of the accretor, and the mass accretion rate, respectively. 
Solving for m, and inserting the measured X-ray luminosity for L x gives, 



where R$ is the radius of the primary in units of 10 9 cm, and R/Rl 1 ~ 0.2 for a system 
like GP Com. For a primary mass, Mi = O.5M ( Morales- Rueda et al. 2003) and R 9 = 1 
(Zapolsky & Salpeter (1969) we obtain rh ~ 8.7 x 10~ 13 M yr _1 . If the accretor is a 
massive O/Ne/Mg white dwarf (suggested by the possibility of oxygen and neon mixing in 
the boundary layer), then the derived accretion rate drops by about a factor of 6. 

This rate is substantially lower than the critical rate above which one expects an optically 
thick boundary layer (Pringle & Savonije 1979; Narayan & Popham 1993). In this regime 
the plasma should have a maximum temperature of « 10 s K, which is similar to the peak 
temperature of 6.5 keV (7.5 x 10 7 K) deduced from our spectral modelling. The presence of 
narrow emission lines also supports spectral formation in a boundary layer. From the width 
of the N VII Lya line, we obtain an upper limit on the velocity dispersion of the emitting 
gas of 250 km s _1 (3a) (This includes bulk and thermal motions). This indicates an origin 
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close to the white dwarf and not, for example, further out in the accretion disk, which would 
have much higher Keplerian velocities. If the white dwarf accretor was synchronized with 
the orbit, then its rotational velocity would be 23 km s" 1 . To have a rotational speed 
equal to our 250 km s _1 limit, a 10 4 km white dwarf would require a spin period of ~ 250 
seconds. Indeed, an upper limit to the Doppler (ie. thermal motion) width for the nitrogen 
line indicates a value close to t250 km s _1 , which supports a slowly rotating, or synchronized 
white dwarf. The narrowness of the X-ray lines also appears consistent with the narrow 
"central spike" component of the helium lines (Marsh 1999; Morales- Rueda et al. 2003), 
which almost certainly originate on the accreting white dwarf. 

The XMM/Newton measurements provide a detailed look at the physics of the boundary 
layer in GP Com. Detailed comparisons with theoretical boundary layer models could provide 
interesting constraints on their temperature, density and rotational profiles. Deeper X- 
ray spectra would provide better temperature and density constraints from the helium-like 
triplets, and likely more lines of less abundant species could be detected. 

The author thanks Richard Mushotzky, Martin Still, Craig Markwardt, and the referee, 
Tom Marsh, for many useful comments and suggestions. 
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Figure 1: A lightcurve of a portion of the XMM/Newton EPIC (PN+MOS) observation of 
GP Com. The binsize is 80 s. Significant variability is evident on a range of timescales. A 
characteristic error bar is also shown (upper left). 
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Figure 2: High resolution spectra of GP Com obtained with the RGS1 (bottom) and RGS2 
(top) spectrometers in the 7 — 32 A band (solid histogram). Identifications are given for 
the prominent lines. The best fitting, variable emission measure plasma model (cevmkl in 
XSPEC) is also shown (solid curve). 
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Figure 3: RGS spectra and best fitting Lorentzian line profiles for the N VI (top) and Ne IX 
(bottom) helium-like triplets. 
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Figure 4: Phase- averaged spectrum of GP Com as observed in the EPIC/MOS detectors 
(MOSl:triangles, MOS2:crosses). Shown are the data and best fitting cevmkl model versus 
wavelength using XSPEC. The fit includes contributions from helium, nitrogen, oxygen, 
neon, sulphur and iron. 
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Table 1: Emission Line Properties for GP Com 



Species 




A (A) 






Flux 






r 


i 


f 




(10~ 4 cm" 2 s- 1 ) 




N VII 




24.781 






1.86 




N VII 




20.910 






0.28 




N VI 


28.792 


29.075 


29.531 


0.34 


0.22 


0.07 1 


Ne X 




12.134 






0.54 




Ne X 




10.239 






0.16 




Ne IX 


13.447 


13.550 


13.697 


0.20 


0.11 


0.09 



1 Not detected, la upper limit. 



